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Time evolution of the probability of quasilocalized massive modes decay trough the tunneling into 
an extra dimension is analyzed in the Randall-Sundrum braneworld model. 
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I. INTRODUCTION 

The idea that our world is a three brane embedded 
in a higher dimensional space-time with non-factorizable 
warped geometry has been much investigated since the 
appearance of papers P, 0. The first solution for a 
warped metric was provided in P]. In this paper it was 
suggested that gravitational interactions between parti- 
cles on a brane in uncompactified five-dimensional space 
could have the correct four-dimensional Newtonian be- 
havior, provided that the bulk cosmological constant is 
fine-tuned to the brane tension. As it was shown in 2] 
the warped metrics could provide an alternative to com- 
pactification of extra dimension used for a possible ex- 
planation of the hierarchy problem 3]. Brane solution 
described in 0, Q has the form 



-2k\z\ 



(1) 



where the parameter k is determined by the bulk cosmo- 
logical and five-dimensional gravitational constants re- 
spectively. An important issue in the study of braneworld 
models is the localization of standard model fields on the 
brane. This problem for the brane |^ was addressed 
in 4J. It was shown that the massless scalar field has 
a zero mode normalizable with the weight exp(— 2fc|z|). 
When this field is given mass the bound state becomes 
metastable against tunneling into an extra dimension . 
Fermion fields are not localized on the positive tension 
brane by gravitational interactions Q- In order to lo- 
calize fermions one can use the method of localization 
described in paper But the adding of bulk mass 
term makes the fermion bound state unstable as well 
trough the barrier penetration into an extra dimension 
0. Based on this phenomenon the possibility of decay- 
ing cold dark matter into an extra space was considered 
in the framework of braneworld cosmology, see @ . The 
comoving density of cold dark matter was taken to decay 
over time with a rate pa? exp(— Ft) corresponding to the 
exponential decay of a metastable state. But in Quan- 
tum Mechanics it is well established that an exponential 
decay cannot last forever if the Hamiltonian is bounded 



below and cannot occur for small times if, besides that, 
the energy expectation value of the initial state is finite 
. From this point of view we consider the dynamics of 
massive matter disappearance on the brane. 



II. SCALAR FIELD 

Due to four-dimensional Poincarc invariance of the 
metric every fields in this background can be decom- 
posed into four-dimensional plane waves 



4> oc exp(ip^a:^)(p(z) , 



(2) 



where coordinate four momentum coincides with the 
physical momentum on the brane. It is natural to assume 
that the three momentum p does not change the picture 
of tunneling into an extra dimension. For the sake of 
simplicity we put p = in what follows. Under this 
assumption the equation 



-2k\z\ ,,2 



M> = , (3) 



gives 



-9,V + 4fcsgn(z)9,<^ + ^2^ = e2'=l"lsV , (4) 

where E = pq. The continuum spectrum of eq. Q starts 
from ^2 ^ 



a{E)J, (f ^'''"') +b{E)Y, 



E 



.k\z 



where v = + 5, 6]. The normalization condi- 

tion 



dze~^^\''\ip{E, z)(f{E', z)^5{E-E') , 



as well as the boundary condition dz^\z=o — can be 
satisfied by taking 
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a{E) 



A{E) 



^l + A^{E) 



, b{E) 



^l + A^(E) 
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A{E) 



Y,^,{E/k) -{u- 2){k/E)Y,{E/k) 



J,-i{E/k) -{v- 2){k/E)J,{E/k) • 
The general solution to Eq.JS)) can be written as ^3 



(j){t, = J Gi{t, z, z')0(O, z')dz' 

— OO 

oo 

+ / G2{t, z, z')(P{0, z')dz' , (5) 



where 



Gi{t, z, z') = e^^''!^'! / dEcos{Et)(p{E, z)ip{E, z) 



G2{t, z, z') = e-^'^t^'l / dE^^^^v{E, zME, z') . 

To compute 4>{t, z) one needs to know the initial data 
(^(0, z), 0(0, z). In what follows we take the following 
condition (/)(0, z) — iEQ(j){0, z) corresponding to the free 
particle on the brane 10) with energy Eq. Correspond- 
ingly the Eq.© takes the form 



z)= / [Gi(<, z, z')+iEoG2{t, z, z')](j){0, z')dz 



(6) 

At i = the particle is confined on the brane. The prob- 
ability of the particle to remain on the brain at instant t 
is given by 

\{m\m? , 

where (/)(0) denotes the brane localized initial state 
and the scalar product is understood with the measure 
exp(-2A:|z|) 



(mm)) = J dze-^'^^^no, z)<l>{t, z) . 

— oo 

From Eq.© one finds 



Re{cl){0, z)\(l){t, z)) = Re 



OO 

J dEe-'^'\G(E)Y 



(7) 



Im{(j){Q, z)\(j){t, z)) = -Eolm 



-lEt 



dE- 



E 



\C{E)\' 



(8) 



ImE 



ReE 




where 



Poles 



FIG. 1: The integration contour. 



In the above E is real variable, but in the following 
it is understood that E can be complex. The function 
|C(i?)| has the simple poles determined by the equation 



i ± A{E) = 



The equation 



A{E) 



k^'^-^k 



+ (2-i.)i/« 



under assumption fj,/k,E/k <^ 1 gives the pole in the 
fourth quadrant of a complex E plane, E — Eq ~- iT, 
where 



Eo 



V2 



Eo 8\k 



Since, when t > and \E\ oo, ^ in the fourth 

quadrant, for evaluating of integrals entering Eqs.([7||Hl) 
one can deform the integration counter as it is shown in 
Fig.l. In this way one finds 

oo 

J dEe~'^'\CiE)\'' = 



oo 

residue term - i J dEe-^* \G{-iE)\^ . (9) 



From this equation one sees that the second term 
gives an imaginary contribution and correspondingly 
Re{4>(0, z)\(j){t, z)) is determined by the residue term 
only. In the same way one finds that Im {4>{0, z) z)) 
is also determined only by the residue term. Following to 
the quantum mechanical formalism, and by taking into 
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account that /^/A:, E^/k <C 1 and the localization width 
of the scalar is ~ k^^ , one can simply take 



In terms of the right- and left-handed components 



MO z)-{ ^/^/l^^ ' for \z\ < (2fc)-i , 

" \ , for \z\ > (2fc)-i . ^^"^ 

Thus, the probability of finding the particle on the brane 
equals one when t = 0. By evaluating the residue term 
and using Ea. l|lU|l one finds 

|((/.(0)|</.(<))|' = 0.04x 



III. SPIN 1/2 FIELD 

The minimal representation of spinors in five dimen- 
sions can be chosen to be four dimensional. The five- 
dimensional Minkowskian gamma matrices can be chosen 
as follows 



= -ir 
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fj'^ 

CT^ 



7 = 



I 

-I / ' 



where cr'' = (1 , ct) , cr^ = (1 , — (?) with ct^ the three 
Pauli matrices. The Dirac equation in the background 

ds^ = e-'^^^^^ri^^dx^'dx" ~ dz^ , 

for the fermion coupled with the domain wall $ reads 



zr^(9, - 2a'{z)) + ie'^f^^r^S^ - g$ 



^ = 0. (11) 



Solutions to this equation are linear combination of wave 
functions of the form 

V' = exp(-ip^a;^)x . 

We do not expect the three momentum p to affect the 
overall picture of particle tunneling into an extra space. 
Therefore, as in the case of scalar field, we assume p = 0. 
Under this assumption the Ea. (|ll|) takes the form 



where 



(ido ~H)ij^O, (12) 



i7 = e-^[7V-7V(9.-2(T')] . (13) 



So that the time evolution of initially brane localized 
fermion ^/'(O, z) is given by 



i;{t, z)^ J dEe-'^\iE, z) {x{EmO)) , (14) 



where the scalar product is understood with the measure 



Hx^Ex 



(15) 



XR 



1 + 7^ 



-X , XL 



1 -7^ 



the Eg. pSfl takes the form 



{d,-2a' -g^)XR = -Ee''XL 



id,-2a' + g^)XL=Ee''XR , 



(16) 



(17) 



where E = pq. After eliminating XR from Eas. (|16l I17f) 
one obtains a second order equation for xl 



+.g$' - {g<^f + E^e^^] xl - . 



(18) 



For the domain wall profile in the thin wall limit one can 
take 

$ = ?;sgn(z) . 

To the right (z > 0) and left (z < 0) of the brane one 
gets 

[dl ~ 5kd, - gkv + 6k^ - g^v^ + E^e^^'-] xl = , 
[dl + 5kd, - gkv + 6k^ - g^v^ + E^e'^'''] xl^Q -(19) 

The solution of Eq. ifT^ reads 



XL 



2k 



a{E)J, ( ^e'^^A +h{E)Y, ('|e'=l^l 



(20) 

where v = {2gv + k)/2k. From Eas. p6ll7|l one finds the 
following boundary condition 

Xl(0) + (5«-2/c)xl(0) = , 

where the prime stands for the derivative with respect 
to |z|. To satisfy this boundary condition as well as the 
normalization condition 

oo 

j dze-^^\^\xL{E, z)xl{E\ z) = 5{E-E') , 

— oo 

the coefficients a{E)^ H-E) should have the form [T]| 

Y,.iiE/k) 



a{E) 



^YS^,{E/k) + Jl^,{E/k) ' 



h{E) = -^1-a^E) 
The solution to Eq.Jinil for E^O 



(21) 
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is localized on the brane as long, as k < 2gv with the 
localization width 2/(2gv — fc) Under assumption 
that the initial state is given by the left-handed particle 
localized on the brane one finds 

oo 

(MOMLit)) = / dEe-^'^'lixLiEMLmi^ ■ (22) 



The function |(xl(£')|V-'l(0))|^ niay have the poles En in 
the fourth quadrant of complex E plane determined by 
the zeros of Hl^}^{E/k)Hl^\{E/k) indicating the pres- 
ence of resonances in the spectrum of KK modes. For 
evaluating the integral H22() one can deform the integra- 
tion contour as it is shown in Fig.l. In this way one 
gets 

(MmLit)) = 

oo 

J2 c„e-^^"* - I f dEe-^'\{xL{-iE)\MO))\' , (23) 

n 

where 

Cn^-27ri lim iE-En)\{xL{E)\MO))\'' ■ 

E — >i?Ti 

So that the probability to find the particle on the brane 
at time t is given by 

|(V'L(0)|V'L(t))|^ - |c„pe2-f'"(-BJt4.interference terms . 

71 

One can easily check that 

d{M0)\Mt)) ^0 

dt t=o 

Namely, by using Eas. H12ll3|l and taking into account 
that 0(0, z) is an even function of z one sees that in this 
expression the odd function is integrated over a symmet- 
ric region. So that the evolution of the decay process 
requires some time to reach the regime of exponential 
decay. For large values of time t ^ k~^ the second term 
in Eq. \2'6\ becomes dominant for which in this limit only 
modes with E k are relevant. Therefore, the non- 
escape probability behaves asymptotically as 

\{MmL{t))? - {kt)-^"^^ . 

Thus, after a long time the decay proceeds with respect 
to the power law. As it is shown in one can choose 
the index v in such a way {2 < v < b) the Hankel function 
H^^\{ijj) to have exactly one zero in the fourth quadrant 
of complex Lo plane. This is the only zero of the product 



h]j'\{uj)H^}_-^{lij) located in the fourth quadrant since 

for real values of Hu^\u!)* = H^\uj*). From the 
trajectories of zeros presented in one sees that if v 

is close to 2.5 {v > 2.5) the function H^\{E/k) has 
the zero E — E^) — iV m the fourth quadrant such that 
Eo/k <^ 1 and T ^ k. Using Ea.f l^ for the initial state 
of the low lying KK resonances one can simply take 

^^(0) ^ / V^^^^''""^^'^' ' 1^1 ^ (2^^^ - ' 

lo, \z\>{2gv^k)-' . 

(24) 

So that at t — the particle is known to be on the brane, 
\z\ < {2gv — k)^^, with probability one. Using Ea. l|24() 
and the previous consideration one can perform the con- 
crete calculations for different low lying KK resonances. 



IV. DISCUSSION AND CONCLUSION 

From the present consideration one sees that if the 
transverse equation contains the second order time 
derivative the corresponding decay law has the expo- 
nential form as long as F ^ i?o. So that under this 
assumption the decay of massive initially brane local- 
ized metastable modes for scalar, vector and gravita- 
tional fields follows to the exponential law. It should 
be stressed that the choice (j){0, z) — iEo<j){Q, z) which 
is essential to this result is not unique. However we find 
this condition to be natural as it corresponds to the free 
particle on the brane |(2Jl with energy Eq corresponding 
to the metastable state. In contrast, the decay of massive 
quasilocalized fermion modes does not follow completely 
to the exponential law. As in the standard quantum me- 
chanical case mentioned in the introduction, initially the 
decay is slower than exponential, then comes the expo- 
nential region and after a long time it obeys a power 
law. The present approach allows one to evaluate the 
preexponential term as well as the whole dynamics of 
the decay and can be important to evaluate dynamics for 
any process based on this phenomenon, for instance to 
examine more precisely the cosmological constraints on 
the paradigm of decaying dark matter in the braneworld 
model. 
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